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ABSTRACT 

The observed relation between supermassive black hole (SMBH) mass (M.) and bulge stellar velocity 
dispersion (cr*) is described by logM. = a + (3 log (cr* /200km s _1 ). As this relation has important 
implications for models of galaxy and SMBH formation and evolution, there continues to be great 
interest in adding to the M. catalog. The "sphere of influence" (r^) argument uses spatial resolution 
to exclude some M, estimates and pre-select additional galaxies for further SMBH studies. This Letter 
quantifies the effects of applying the argument to a population of galaxies and SMBHs that do not 
follow the M, — cr* relation. All galaxies with known values of cr*, closer than 100 Mpc, are given 
a random M, and selected when n is spatially resolved. These random SMBHs produce a M. — cr* 
relation of a = 8.3 ±0.2, f3 = 4.0 ±0.3, consistent with observed values. Consequently, future proposed 
M. estimates should not be justified solely on the basis of resolving r%. This Letter shows the observed 
M. — cr* relation may simply be a result of available spatial resolution. However, it also implies the 
observed M, — cr* relation defines an upper limit. This potentially provides valuable new insight into 
the processes of galaxy and SMBH formation and evolution. 

Subject headings: black hole physics — galaxies: bulges — galaxies: fundamental parameters 



1. INTRODUCTION 

There has been a growing database of direct super- 
massive black hole (SMBH) mass (M.) estim ates from 
the ce nters of nearby galactic bulges (e.g., iGrahaml 
2008b). While the limits of our current abilities to sig- 
nificantly expand this database may have been reached 
([Batcheldor fc Koekemoeri [2009), the last decade has 
seen a wealth of M. estimates t hat have increased 
the SMBH catalog from 13 or 26 (iFerrarese fc Merrittl 
20001: IGebhardt et all 120001) t o ~70 (jGraharX l2008b: 
Hul 120081: iGultekin et all 120091) . An intense interest in 



populating the SMBH database was sparked by observed 
correlations between M, and fundamental properties 
of t hei r host bulges (e .g.. i Kormendv fe Richstonc 
19951 : iMagorrian et al.l Il998t iFerrarese fc Merritt 



2000[ IGebhardt et all 



Graham et al.l 



Ferrarcsc 2002; Ma rconi fc Hunt 2003: B aes et al.l[200l 



Haxing fc Rix 2004; Pizzclla et al. 2005). These M. seal 



ing relations h ave generated numerous theoretical inves- 
tigat i ons (e.g., Ciotti fc van Alb ada 2001; Adam s~et al.l 
120031: ICattaneo et~aT1 120051: iRobertson et all 120061 ) and 
have possibly added valuable limits to evolutionary 
models (e.g., iHeckman et all 120041 : IWvithe fc Loebl 
120051: iTreu etalll2007l : ICiotti 200l . 

The degree to which a SMBH's sphere of influence, 
ri, is resolved has been used a s a quality measure for 
M, estimates (jF^ rrareseJ 120021 : iMarconi fc Huntl 120031 : 
IValluri et al.ll20 04). The only method available to a pri- 
ori determine i f a M, estima te can be made is to assume 
ri = GM,/a1 ([Peeb les 1972). However, to calculate in 
galaxies with known cr* , M, is estimated using the M m — 
er» relation given by logM. = a + /31og (cr* /200km s _1 ), 
where a = 8.1- 8.2 and /? = 3.7-4 9. The observed scat - 
ter, e, is 0.4 dex ([Novak et al.ll2006l : iGraham fc Lill2009f ). 
Following this, Figure [T] demonstrates where will be 
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resolved, given a spatial resolution of Sft = O'.'l. A value 
of K = O'.'l is used here as that is the typical FWHM 
of the HST PSF. To date HST has been responsible for 
most M, estimates. 

IFerrarese fc Fordl ([2005L FF05) discussed the l i mited 
abilities of HST to resolve r^, and IGultekin et all (|2009l 
G09) found the M. — cr* relation to be biased when ap- 
plying the ri argument. The influence of rj cuts on the 
M % — cr* relation is continued in this Letter with two sig- 
nificant advances. First, a sample of all galaxies with a 
known cr* (<100 Mpc) is used. It is important to only 
use these galaxies as there are no cr*=400km s _1 galaxies 
at 1 Mpc, for example. Second, random M, estimates 
are applied to each galaxy, i.e., no galaxy is assumed to 
intrinsically fall on the M, — cr* relation. This ensures no 
pre-selection of galaxies that lie on the M. — cr* relation. 
Throughout, a distinction is made between the observed 
M, — cr* relation (published values) and the observable 
M. — cr* relation (the relation that can be fitted using 
the data simulated here). 

2. METHODS AND RESULTS 

In short, these simulations take a galaxy with a known 
cr* and distance, assign a M,, then calculate r^. The 
argument is then applied; if r t is unresolved the galaxy is 
removed from the sample (the low mass cut). If the as- 
signed M, generates a galaxy that lies above the observed 
M, — cr* relation, it is also removed from the sample (the 
high mass cut). The observable M,—o~* relation is then fit 
to the remaining galaxies using a Levenberg-Marquardt 
least-squares add-on to IDL. 

An SQ L search of the H yperleda 1 catalog was per- 
formed ([Paturel et al.l [2003). All galaxies with a kine- 
matical distance modulus less than 35.0 (100 Mpc) were 
found. Of the 49740 galaxies returned, 2518 have pub- 
lished values of cr*. The incompleteness of this sample, 

1 http://leda.univ-lyonl.fr/ 
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a, (km s ■') 

Fig. 1. — Observable M m — ct„ plane at discrete distances based 
on 3? = 071 and a = 8.2,0 = 4.6 (solid line). In areas above the 
dashed distance lines, r% will be resolved if the galaxy with the 
given a, hosts the given SMBH. 

due to difficulties in measuring c* in faint galaxies, is 
noted. Values of (in arcsecs) were then calculated for 
each galaxy by assigning a M.. 

In the first cases, the low mass cut was made by apply- 
ing 5R = O'/l to represent the HST PSF. However, several 
values of 5ft are applied later. The high mass cuts were 
then made by considering the observed distribution of 
galaxies in the M, — cr* plane; there is a clear absence 
of over-massive SMBHs, i.e., there are no <7*=50km s _1 
dwarf spheroids in the Local Group with M.~ 10 9 Af Q . 
Such a nearby over-massive SMBH population would 
have been detected, and therefore it is assumed that these 
objects do not exist. Galaxies selected for the high mass 
cuts were determined by assigning a scatter (e = 0.4 dex) 
and upper limits to the M. — cr* relation, (M. — er») u . 
Galaxies with ri > r max (M max ) were removed from the 
sample, where M max = [e + 10 Q " (cr/200km s- x )^]M . 

Each galaxy was first assigned 90 separate values of 
M. using a step size of 10 01 from 10 10 Af Q to 10 10 M Q . 
In this case, the high mass cut was made by apply- 
ing the a u — 8.1, (3 U = 4.2 relation of G09 and the 
a u = 8.2, /3 U = 4.9 relation of FF05. Figure [2] shows 
the distribution of observable galaxies in the M. — <r» 
plane assuming every galaxy within 100 Mpc can host 
any M. and that 5ft = O'/l. The fit to this sample (red 
line) is a = 8.4, /3 = 3.5. 

However, galaxies host single (or binary) SMBHs, 
therefore it is shrewd to assign each galaxy a single 
(uniformly sampled) random value for M, in the range 
10 1 — 10 10 M Q . The distribution of observable galaxies 
in the M m — <t* plane using a random sample of M. 
is shown in Figure [3] using the red open circles. In 
this case, the low cut was made using 5ft = O'/l, and 
the high mass cut was made using a (M, — cr»)„ re- 
lation of a u = 8.1, P u = 4.2, e = 0.4 (G09). The 
M. — (7* relation fit to these observable galaxies (red line) 
is a u = 8.3, j3 u — 4.1, e = 0.2 dex. An observable M. — a* 
relation of a u = 8.3, j3 u — 4.6, e = 0.4 dex is found using 
the high mass (M. — er*) u relation of FF05. Figure |S] 
also shows all observed gala xi es based o n the combined 
catalogs of iGrahaml (|2008bf ). |5ul <|2008[) and G09. No 
distinction is made between "good" and "bad" M, es- 
timates, or differences in quoted <r*; all estimates are 
plotted (143 total). 

The similarity between the M m — a* distribution of 
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Fig. 2. — Simulated M. — <r» data points (black circles) based 
on a SR = O'/l criteria. This demonstrates the observable region of 
the Mm — a* plane with the assumption that there is no M, — ct„ 
relation. The fitted M, — <r, relation (red line) is a = 8.4, = 3.5. 



TABLE 1 

Effects of the Cutoff From Averaged 
Random M. Distributions 





Observed ct 


Observed 


Observed e 




(Zero-point) 


(Slope) 


(Scatter) 


Variable Upper Cutoffs 


0705 


8.25 ± 0.17 


3.80 ±0.33 


0.41 dex 


O'.'IO 


8.25 ±0.21 


3.90 ± 0.34 


0.27 dex 


0715 


8.25 ±0.23 


3.93 ± 0.34 


0.22 dex 


0720 


8.25 ± 0.27 


4.01 ± 0.38 


0.18 dex 


Fixed Upper Cutoff 


0705 


8.44 ± 0.03 


3.93 ±0.23 


1.82 dex 


0710 


8.57 ±0.02 


4.02 ±0.21 


1.22 dex 


0715 


8.64 ± 0.02 


4.03 ± 0.25 


0.97 dex 


0720 


8.68 ± 0.03 


4.07 ±0.32 


0.85 dex 



Note. — Average results from applying vari- 
able and fixed M. — it* upper limits to 50 random 
M, samples. A fixed upper limit of a u = 8.7 and 
P u = 5.0 was used based on the upper limit fit in 
Figure [4] If there was no intrinsic M m — <t„ rela- 
tion, these are the values of the M. — ct* relation 
that would be observed using the different spatial 
resolutions listed. 



observable random mass SMBHs and observed SMBHs 
masses is striking. However, as this result derives from 
a single random sampling of M # , it represents a single 
possible observable M m — er* relation. If galaxies intrinsi- 
cally have random M % values, then the range of observ- 
able SMBHs could have a distribution given by Figure [5J 
Therefore, 50 separate random M» samplings were then 
made for each galaxy, and in each case a fit to the ob- 
servable M, — a* relation was performed. In addition, 
as the previously used value of — O'/l only applies to 
HST observations, the analysis was repeated for a range 
of Ti lower cut offs. Finally, as the high mass cuts may 
not actually be defined by the G09 and FF05 fits, values 
of a u from 7.8 to 8.7 (in 0.1 steps) and /3 U from 3.6 to 
5.4 (in 0.2 steps) were used to create an addition 100 
individual (Aim — cutoffs. The mean values of a, /3 
and e, derived by applying these different 5ft criteria, and 
from using these variable upper limits, are presented in 
Table [TJ 
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a, (km s ■') 

Fig. 3. — Distribution of observable galaxies from a random 
distribution of M, (red circles) using the (M. — <7*) u relation of 
G09. All present estimates of Mm and tr* are shown as filled black 
circles, with uncertainties. The fit to the observable galaxies (red 
circles and line) is a = 8.3, (5 = 4.1 
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a, (km s ■') 

Fig. 4. — Fitting the observed (M. — cr*) u relation. Filled black 
circles are the M, — tr* points, those with error bars are the points 
used for the upper fit (shown in red). The red upper limit fit is 
a = 8.7, ft = 5.0. The two solid black lines are the observed G09 
and FF05 fits. 



Table Q] has several notable features. First, as 5ft rises 
e falls. This is expected because the range of observ- 
able SMBHs decreases with larger values of 5ft. Sec- 
ond, the slope (/?) remains consistent. Finally, in the 
case of the variable upper limits, the zero-point (a) re- 
mains constant and is consistent with the mid-value of 
the (M. — cr*) u limits imposed on the sample. To test 
whether these values of a are a consequence of the high 
mass cut conditions, the analysis was repeated using a u 
limits of 7.0 and 8.8 (in 0.2 steps). Applying a 5ft = O'.'l 
low mass cut, an observable relation with a — 8.0 ± 0.5 
was found. Again, the values of a are consistent with the 
mid-value of the imposed limits. However, an estimate 
of the observed (M. — o~*) u relation can be made. The 
observed SMBH sample in Figure|3jwas used to define an 
upper limit by fitting all observed SMBHs that fall above 
the FF05 and G09 relations, and above e = 0.4. The fit 
to these over-massive SMBHs produce an upper limit of 
a u = 8.7, p u = 5.0 (Fig. gj. This (M. - tr*) u fit was 
then applied to 50 random values of M, in each galaxy, 
and the observable M, — tr* relations were then fitted in 
each case. The average of these fits are also presented 
in Table [TJ In this case, the fitted values of a and f3 are 
consistent with the variable upper limits, but the values 
of e are significantly larger. This is expected, as these 
fits were applied to data that had the maximum range 
in the M, — tr* plane due to the high observed values of 
a u and j3 u . 

3. DISCUSSIONS 

In summary, as a result of applying the r*j argu- 
ment, a M, — tr* relation consistent with observed values 
(a w 8, f3 w 4, e « 0.3 dex) can be fitted to a sample 
of galaxies that contain random mass SMBHs, and as a 
consequence do not follow a M m — tr* relation. The r, 
argument removes low mass SMBHs where r, would not 
be resolved, and high mass SMBHs where Ti should have 
been resolved if such a population were present. There- 
fore, for scaling relations to be of value in constraining 
galaxy evolution models, M % estimates must not be solely 
proposed on the basis of critically resolving values of r, 
derived from cr*. 

It is unlikely that an intrinsic M, — cr* relation has 
been observed as a result of critically sampling in the 



M m — tr* plane, and it will be some time until we are able 
to sample a significant area rightward of the observed 
Mm — tr* relation (e.g., FF05: iBatcheldor fc Koekemoerl 
2009). Therefore, it is clear that caution must be used 
in the application of the observed M, ~ o* relation until 
this issue can be resolved. 

The possibility of there being no tight intrinsic M % — tr* 
relation, and that the distribution of observed galaxies 
in the M. — cr* plane determine an upper limit only, is 
now discussed. As the same M, estimates are used to 
define the other M, scaling relations, this also implies 
that all observed scaling relations are upper limits. There 
are a number of questions arising from this possibility, 
beginning with the most fundamental. Why do observed 
SMBHs fall on the M, — cr* relation? There are several 
SMBHs with such high quality M, estimates (Sgr A*, 
NGC 4258, M87), that there clearly is a relation between 
Mm and tr* in some galaxies. However, the M. — cr* 
plane has been increasingly populated with less certain 
Mm estimates that have potentially been included on the 
assumption that ri is spatially resolved, i.e., that they 
follow a previously estimated M, — cr* relation defined 
by higher quality M, estimates. In these cases, as the 
simulations presented here show, an observed M m — cr* 
relation will arise simply as a result of the selection 
effect, even if M. is randomly distributed within galaxies. 

Does a population of under-massive SMBHs exist, 
and have they been detected? If the M. — tr* rela- 
tion is an upper limit, then there should be galaxies 
that host low mass SMBHs, as sugge sted by simulations 
(jVittorini et al.ll2005HVolonteri| [2007'). Indeed, both the 
simulations presented here, and the observations plot- 
ted in Figure [31 show that under-massive SMBHs can 
be, and have been, detected. Some not able cases are 
that of NGC 4 435 (ICoccato et alJ l20M ). a sample of 
barred gal axies (lGrahamll2008al) a nd narrow-line Seyfert 
1 galaxies (M athur fc Grupe 2005). In fact, the evidence 
to suggest the presence of under-massive SMBHs is not 
matched by any evidence for a significant population of 
over-massive SMBHs leftward of the observed M, — cr* 
relation. 

It is important to note some intricacies with the two 
dominant techniques for measuring M. (stellar and gas 
dynamics). In the case of stellar dynamics, there are 
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systematics to the mod els that may allow a large range 
of M. in a given bulge (|Valluri et al.l 120041) . In the case 
of gas dynamics, it is unclear w hat the inclination o f 
the nuclear gas disk may be (e.g., iMarconi et al.l 120031) . 
Both methods allow the potential for many of the current 
M. estimates to in fact be upper limits, i.e., the true 
M. — cr* plane may have a large distribution of under- 
massive SMBHs. Therefore, under-massive SMBHs may 
have been observed, their mass over-estimated, and their 
impact over-looked. 

All SMBH mod els allow stringent upper limits to b e 
placed on M. (e.g. JSarzi et al.ll2002HBeifiori et al.ll2009f ). 
However, these M. upper limits will also be dependent 
on the available spatial resolution. Any kinematical data 
will have two velocity points spatially separated on a 
scale of 5i. An upper limit to M, is estimated by in- 
cluding an increasing dark mass until the derived model 
becomes inconsistent with the data, i.e., a higher upper 
limit to M.will be estimated using a lower 5ft. There 
are still important constraints that can be added to the 
M. — cr* plane from estimates of M m upper limits, how- 
ever, as upper limits that fall below the observed M m — cr* 
relation provide the same evidence for under-massive 
SMBHs as would a tightly constrained low mass M,. At 
present, most M. upper limits are based on data derived 
from gas dynamics. Consequently, these limits generally 
fall above the observed M, — er* relation due to unknown 
amounts of line broadening from non-gravitational pro- 
cesses, and uncertainties in the inclination if the nuclear 
gas disk. 

If the M, — cr* relation represents an upper limit in 
the M. — cr* plane, then what is this limit? In § [2] an 
upper limit of a u = 8.7, f3 u = 5.0 was found based on 
the distribution of observed SMBHs leftward of the ob- 
served M, — cr* relation. This is likely a good approxi- 
mation to the (M, — cr*) u limit as there are no reports 
of a steeper relation. However, this estimate does not 
include the potential ov er-massive SMBHs fr om the up- 
per limits calculated by iBeifiori et al.l (|2009f ). Including 
these limits to the upper limit sample gives values of 
a u = 8.8, (3 U = 3.9 to (Al. — a*) u . As expected, due to 
the addition of SMBH limits at lower M,, this (M, — cr*) u 
limit is more shallow with a higher zero-point. Including 
these limits at the lower M. end of the M, — cr* plane 
addresses, in part, a limitation of the sample used here. 
As already noted, the cr* catalog used here is likely in- 
complete due to the difficulty of measuring cr* in faint 



galaxies at greater distances. In addition, the cr* catalog 
likely contains inhomogenuous measurements that may 
not translate from bulge to bulge. 

What are the consequences to galaxy evolution models 
if there is only a (M, — cr*) u relation? First, galaxies will 
no longer be required to obey the M. — cr* relation, and 
could host a SMBH with any M. below (M.-cr*)„. Mod- 
els that include feedback from the SMBH to the galaxy 
will then need to be carefully reconsidered. While the 
SMBH will undoubtedly have some influence on a por- 
tion of the host galaxy, it would not need to affect large 
scale properties; evolution of the SMBH would be a result 
of host galaxy evolution. An upper limit in the M, — cr* 
plane would also represent the pinnacle of SMBH evolu- 
tion as a function of cr* , in which galaxies evolve up to the 
(M. — cr*) u limit. A signature of such a scenario could 
be a cosmic variation in e (the scatter would increase 
with redshift) and an observed M m — cr* relation that 
does not exceed (M, — cr*) u . If the distribution of M, is 
random within bulges, then when compared with the lo- 
cal observed M. — cr* relation, M, estimates from higher 
redshi ft could fall to the left or the right. iTreu et all 
(2007) find a population of z=0.36 Seyfert 1 galaxies 
that lie above the local observed M, — cr* relation by 
A log cr = 0.13, AlogM.= 0.54, but this population still 
lies below the (M, — cr*) u limit estimated here. Finally, 
if SMBHs can reside anywhere below the (M. — ct») u 
limit, then the local black hole mass function may have 
been over-estimated. This would relax the observation 
that merging is not important and that SMBH gro wth is 
dominated by accretion (e.g.. IMarconi et"aL| [2004V This 
potentially allows anti-hierarchical SMBH growth to no 
longer present problems for hierarchical galaxy formation 
models. 
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